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An implementation method ol a gate in a quantum computer is studied in terms of a finite number 
of steps evolving in time according to a finite number of basic Hamiltonians, which are controlled 
by on-off switches. As a working example, the case of a particular implementation of the two qubit 
computer employing a simple system of two coupled Josephson junctions is considered. The binary 
values of the switches together with the time durations of the steps constitute the quantum machine 
language of the system. 
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I. INTRODUCTION 



In classical computing the programming is based on 
commands written in the machine language. Each com- 
mand is translated into manipulations of the considered 
device, obtained by electronic switches. In quantum com- 
putation quantum mechanics is employed to process in- 
formation. Therefore the conception of the quantum 
computer programming, the structure of the quantum 
machine language and the commands are expected to be 
quite different from the classical case. Although there 
are differences between classical and quantum comput- 
ers, the programming in both cases should be based on 
commands, and a part of these commands is realized by 
quantum gates. Initially, one of the leading ideas in quan- 
tum computation was the introduction of the notion of 
the universal gate Q . Given the notion of a universal set 
of elementary gates, various physical implementations of 
a quantum computer have been proposed Naturally, 
in order for an implementation to qualify as a valid quan- 
tum computer, all the set of elementary gates have to be 
implemented by the proposed system. This property is 
related to the problem of controllability of the quantum 
computer. The controllability of quantum systems is an 
open problem under investigation |^ . Recently, |^ at- 
tention was focused on the notion of encoded universality, 
which is a different functional approach to the quantum 
computation. Instead of forcing a physical system to act 
as a predetermined set of universal gates, which will be 
connected by quantum connections, the focus of research 
is proposed to be shifted to the study of the intrinsic 
ability of a given physical system, to act as a quantum 
computer using only its natural available interactions. 
Therefore the quantum computers are rather a collection 
of interacting cells (e.g. quantum dots, nuclear spins, 
Josephson junctions etc). These cells are controlled by 
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external classical switches and they evolve in time by 
modifying the switches. The quantum algorithms are 
translated into time manipulations of the external clas- 
sical switches which control the system. This kind of 
quantum computer does not have connections, which is 
the difficult part of a physical implementation. Any de- 
vice operating by external classical switches has an in- 
ternal range of capabilities, i.e. it can manipulate the 
quantum information encoded in a subspace of the full 
system of Hilbert space. This capability is called en- 
coded universality of the system [Q . The notion of the 
encoded universality is identical to the notion of the con- 
trollability 1^ of the considered quantum device. The 
controllability on Lie groups from a mathematical point 
of view was studied in |, |, g, [o|. The controllability 
of atomic and molecular systems was studied by several 
authors see review article |^ and the special issue of the 
Chemical Physics vol. 267, devoted to this problem. In 
the case of laser systems the question of controllability 
was studied by several authors and recently in ||l^, pT| . 
In the present paper we investigate the conditions in or- 
der to obtain the full system of Hilbert space by a finite 
number of choices of the values of the classical switches. 
As a working example we use the Josephson junction de- 
vices in their simplest form jl^, |l^, |l^ , but this 
study can be extended in the case of quantum dots or 
NMR devices. 

In this paper the intrinsic interaction of a system op- 
erating as a universal computer is employed instead of 
forcing the system to enact a predetermined set of uni- 
versal gates 0. As a basic building block we use a sys- 
tem of two identical Josephson junctions coupled by a 
mutual inductor. The values of the classical control pa- 
rameters (the charging energy Ec and the inductor energy 
El) are chosen in such a way that four basic Hamilto- 
nians, Hi,{i = 1, ... ,4) are created by switching on and 
off the bias voltages and the inductor, where the tunnel- 
ing amplitude Ej is assumed to be fixed. Our procedure 
allows the construction of any one-qubit and two-qubit 
gate, through a finite number of steps evolving in time 
according to the four basic Hamiltonians. Using the two 
Josephson junctions network a construction scheme of 
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four steps for the simulation is presented in the case of 
an arbitrary one-qubit gate, and of fifteen steps, for the 
simulation of an arbitrary two-qubit gate. The main idea 
of this paper is the use of a small number of Hamiltonian 
states in order to obtain the gates. 

Our results can be formulated in terms of a quantum 
machine language. For example, each two-qubit gate U 
is represented by a command of the language containing 
15 letters i.e. the fifteen steps mentioned above. Each 
letter consists of a binary part (the states of the on-off 
switches), which determines the basic Hamiltonian used, 
and a numerical part corresponding to the time interval. 

Our proposal can be generalized for 7V-qubit gates 
{N > 2) belonging to the SU{2^). In that case iV + 2 
basic Hamiltonians are needed to implement such a gate. 
The generalization of these ideas is under investigation. 

The paper is organized as follows: In section II we 
present for clarity reasons the formalism of Josephson 
junctions one-qubit devices. In section HI we apply the 
formalism to two qubit gates and simulate two-qubit 
gates and the possible 7V-qubit generalizations are dis- 
cussed. In sectiori^^ the results are summarized. Finally 
in the Appendix Hthe essential mathematical feedback 
is provided. 



II. ONE-QUBIT DEVICES 



1, while all other coherent charge states, having much 
higher energy, can be ignored. These charge states cor- 
respond to the spin basis states: 
I t> corresponding to Cooper-pair charges 
on the island, and 

I |> corresponding to 1 Cooper-pair charges. 
In this case the superconducting charge box reduces to a 
two-state quantum system, qubit, with Hamiltonian (in 
spin 1/2 notation): 



H=^E,a,-^Eja,, a,| t>= | T> 
<Jz\ i>= -I l> ■ 



and 



(1) 



In this Hamiltonian there are two parameters the bias en- 
ergy Ec and the tunneling amplitude Ej. The bias energy 
Ec is controlled by the gate voltage Vg of Fig |l|, while the 
tunneling amplitude Ej here is assumed to be constant 
i.e. it is a constant system parameter. The tunneling 
amplitude can be controlled in the case of the tunable 
effective Josephson junction, where the single Josephson 
junction is replaced by a flux-threaded SQUID ||l^, but 
this device is more complicated than the one considered 
in this paper. 

The Hamiltonian is written as: 



H ^ -AE{rj){cosrja, 



sin TycTj:) 



(2) 



where rj is the mixing angle 



rj = tan 



-I El 

Er^ 



The energy eigenvalues are 



E+^± 



AE{r^) 



FIG. 1: One qubit device 

The simplest Josephson junction one qubit device is 
shown in Fig |^. In this section we give a summary of the 
considered device. The detailed description and the com- 
plete list of references can be found in the detailed review 
paper |jl^, section II] . The device consists of a small su- 
perconducting island ("box"), with n excess Cooper pair 
charges connected by a tunnel junction with capacitance 
Cj and Josephson coupling energy Ej to a superconduct- 
ing electrode. A control gate voltage Vg (ideal voltage 
source) is coupled to the system via a gate capacitor Cg . 

The chosen material is such that the superconducting 
energy gap is the largest energy in the problem, larger 
even than the single-electron charging energy. In this 
case quasi-particle tunneling is suppressed at low temper- 
atures, and a situation can be reached where no quasi- 
particle excitation is found on the island. Under special 
condition described in ref only Cooper pairs tunnel 
coherently in the superconducting junction. 

The voltage Vg is constrained in a range interval where 
the number of Cooper pairs takes the values and 



and the splitting between the eigenstates is: 



AEir^) 



The eigenstates provided by the Hamiltonian (|l|) 
denoted in the following as |-|- > and |— >: 

1+ > = cosf I i> +sin§ I T> 
I- > ^sin§U> -cosfl T> 



are 



(3) 



To avoid confusion we introduce a second set of Pauli 
matrices p = {px, Py, Pz), which operate in the basis |+ >, 
|— >, while reserving the a operators for the basis of | t> 
and I i>: 



pz = 1+ >< +1 
Px = 1+ >< -| 

Py^l\- >< +1 



+ I- X +1, 
- i\+ >< -I 



In the proposed model we assume that the device of Fig 
|l| has a switch taking two values 1 and 0, corresponding 
to the switch states ON and OFF. This switch controls 
the gate voltage Vg, which takes only two values either 
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Vid or Vdog, where the first one corresponds to the idle 
Hamiltonian, while the second one corresponds to the 
degenerate Hamiltonian. 

The idle point can be achieved for a characteristic 
value of the control gate voltage Vg = Vid, corresponding 
to a special value of the bias energy and to the phase 
parameter rj = rji^. At this point the energy splitting 
AE{ri) achieves its maximum value, which is denoted by 
AE. For simplicity reasons we reserve the symbol Ec 
for the bias energy corresponding to the idle point and 
by definition, the Hamiltonian at the idle point then be- 
comes: 



El 



Hk] = -I^CTz —(7x 



\aep. 



(4) 



At the degeneracy point rj — ^ the energy splitting 
reduces to Ej, which is the minimal energy splitting. 
This point is characteristic for the material of the Joseph- 
son junction and corresponds to a special characteristic 
choice of the control gate voltage Vg = Vdog- 



H, 



El 



E 



dcg 



- (sin 77idPz - cos rjidPx) (5) 



The system is switched in the state OFF (or 0) cor- 
responding to the degenerate Hamiltonian (||) during a 
time interval ti , then the system is switched to the state 
ON (or 1) i.e. the idle Hamiltonian (^) during a time 
interval t2 and it comes back to the initial degenerate 
Hamiltonian during the time t^ . The general form of the 
evolution operator is: 



U : 



(6) 



The operators Hi^ and .ffdcg and their commutator 



[Hid , 



ErEj 



dogj 



form a (non-orthogonal) basis of the algebra su{2). 
Therefore the pair Hid, Hdcg generates su{2) by taking 
these elements and all their possible commutators and 
their linear combinations. That means that the combi- 
nation of three terms as in equation (^) for all the triples 
{tiTh^h} cover all the matrices belonging in SU{2). 
Thus we conclude that every 2x2 matrix U in SU{2) 
can be achieved by a device as in Fig ^, with manipula- 
tion of the binary switch permitting to the Hamiltonian 
two possible states i.e. the idle one and the degenerate 
one. The above described manipulations can be codified 
by a rudimentary Quantum Machine Language (QML) 
for the one qubit device. In that elementary language 
the gate U corresponds to a command of the language, 
each command is constituted by (three) letters, each of 
them having the form of a pair 

{e, t}, e = (OFF), or 1 (ON), and < t < oo 



i.e. the command corresponding to equation 
lyzed in the following (at most three) letters: 



U 



{0,ti} 



{1,^2} 



{0,t3} 



The one qubit gates are 2x2 unitary matrices belonging 
to the group U(2). Each element in the group U{2) can 
be projected up to one multiplication constant to an ele- 
ment of group SU{2). Evidently the elements generated 
by the evolution operator (||) belong to SU{2). Through- 
out this paper we shall use projections of U{2^) matrices 
in SU{2^), using the symbol "i— >" to denote this projec- 
tion. Let us consider the fundamental one qubit gates 
or commands NOT, vNOT, Hadamard, Phase Shift and 
their SU{2) projections: 



NOT 



1 

1 



g '4 g'4 



e SU{2), 



\/not = 

\/2 \ e'-4 e '4 
^(1 + ^^.) 



Had 



72 I 1 -1 



V2 



where 



,h _ ,h 

^1 — ''3 



2 ( cos 



Ej + E, 
2 Be 



^sin ^ 



Ej 

AE 



V2Ea 



AE 



PhS 



1 

c^'t' 



where 



,ph 



,ph 



ph 



2 sin" 



Ej 

' AE sin ■ 



Ea 



AE 



The above formulas imply the following analysis of 
these commands in letters (Table ^ . 

The above analysis of a quantum gate in letters is 
rather trivial in the one qubit case and it was presented 
for clarity reasons, but the similar construction is far from 
evident and quite complicated in the A^-qubit case. 
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NOT 



VNOT 



{0,V2igj} 



Had 



PhS 



{Mfy 



TABLE I: Letter Analysis of One-qubit gates 



III. TWO-QUBIT DEVICES 

In order to perform one and two qubit quantum gate 
manipulations in the same device, we need to couple pairs 
of qubits together and to control the interaction between 
them. For this purpose identical Josephson junctions are 
coupled by one mutual inductor L as shown in Fig |^. 
The physics and the detailed description of the coupled 
Josephson junctions are discussed and reviewed in pq . 
For L = the system reduces to a series of uncoupled, 
single qubits, while for L — > oo they are coupled strongly. 
The ideal system would be one where the coupling be- 
tween different qubits could be switched in the state ON 
(or state '1') by applying an induction via a constant 
value inductor L and in the state OFF (or state '0') cor- 
responding to L = and leaving the qubits uncoupled in 
the idle state. The Hamiltonian for a general two-qubit 



to Ec (or ON or '1' state) corresponding to the one qubit 
idle state. Also the parameter takes two values. The 
one is El — {OFF or '0' state) corresponding to a an 
uncoupled two qubit state and the other one has a fixed 
value {ON or '1' state). For the sake of simplicity we use 
the symbol E^ for this induction amplitude. Using this 
combination of parameter values or binary switches val- 
ues, we can obtain the following four fundamental states 
of the Hamiltonian (^: 

Hi: where both of the junctions are in the idle state 



0). 

Hi : 



i?c2 = Ec), while they are uncoupled {El 



i?.(aW+af)-i£;,(aW+4^)) (8) 



This Hamiltonian corresponds to the switches 
choice {Ec,, Ec,, El) ^ (1,1,0). 

H2: where both of the junctions are in the degenerate 



state {Ec, 
coupled. 



Ec, = 0), while the two qubits are 



H, = -lEj{aP+a(^^)-lELai'^ai'^ 
corresponding to the switch choice (0, 0, 1) 



(9) 



E- 



H3: where the first junction is in degeneracy {Ec, ~ 0), 
the second is in the idle state {Ec, — Ec) and they 
are uncoupled {El = 0). 



Hs = lEca^^^ -^Ej (aW 



(10) 



FIG. 2: Two qubit device 



corresponding to the switch choice (0, 1,0). 

Hi. where the first junction is in the idle state {Ec, — 
Ec), the second is in the degeneracy {Ec, = 0) and 
they are uncoupled {El = 0). 



system is written: 

H = Hi + H2 + Hint 



(1) 



hE, ai'^ ai'^ 



(7) 



For an explanation of the formalism used in this section 
see Appendix |^. In the case of two identical junctions we 
have Ej-^ = Ej, =^ Ej, since the tunneling amplitude of 
the junction is a system parameter, depending on the ma- 
terial. Under these conditions the two coupled Josephson 
junctions Hamiltonian will be controlled by the following 
control parameters '-Ec,, Ec,, El, which will be called 
switches. The first two parameters are controlled by the 
gate voltages Vg^ , Vg, , while the last parameter is related 
to the inductor switch L. In the proposed model each of 
the parameters Ec,, Ec, can have two values or Ec- 
The first is the state '0' (or OFF) corresponding to the 
degenerate one qubit state, while the other one is equal 



H. 



lEcai'^ 



1 



Ej {ai- 



(1) 



2 ^ 2 
corresponding to the switch choice (1,0,0). 



(11) 



These four Hamiltonian forms are linearly independent 
and they can generate the su(4) algebra by repeated com- 
mutations and linear combinations. For a detailed discus- 
sion see the discussion in Appendix |^, equation ( [A1C| ). 
From Proposition A. 2, any elementary two-qubit gate, 
which is represented by a unitary 4x4 matrix U gSU(4), 
can be constructed as follows: 



iHsts ^-iH,t, ^-iH,t, 



(12) 



From a physical point of view any two qubit quantum 
gate can be obtained by 15 time steps. At the k-th 
step the device is put at one of the Hamiltonian states 
Hi, H2, H3 or H4 by appropriate manipulations of the 
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switches during a time interval tk, k = 1,2, ...,15. Us- 
ing this fact a Quantum Machine Language (QML) can 
be defined as in section II. Each gate, corresponding 
to the 4x4 matrix U, is associated to a command of 
the QML. Each command consists of at most 15 letters 
or steps, each of them being a coUection of 4 numbers 
{ei, 62, i, t} of the following form: 



letter <^ {ei, 62, t} 



<t<oo 



ei 





if V,, = 


Vdcg = 




= 




= 1 


if Vg, = 


= 


=> E,, 




£2 


= 


if Vg, = 


Klcg = 


^ Ec, 


= 




= 1 


if V,, = 


V^d = 


^ Ec., 


= E, 


i 


= 


if L = 




^ El 


= 




= 1 


if L 7^0 




^ El 


^0 



The gate U given by equation (|T^) corresponds to one 
command, which contains 15 letters and is presented in 
Table Each letter is the codified command which in- 
dicates the state of the binary switches and the time in- 
terval. 



2. The SWAP gate, which interchanges the input 
qubits: 



SWAP 



/ 1 

10 

10 

V 1 . 

/ (1) ( 



/1)^(2) 



- ai^ Vi^^ + I I) 



3. The QFT4 gate, the gate of the Quantum Fourier 
Transform (the quantum version of the Discrete 
Fourier Transform), for 2 qubits. A very useful 
gate for the implementation of several quantum al- 
gorithms (e.g. Shor's algorithm): 



QFT4= \ 




U 



{l,l,0,ti} 



{o,o,i,fe} 



{0,l,0,t3} 



{l,0,0,t4} 



{i,i,o,fe} 



{0,0,1,^6} 



{0,l,0,t7} 



{1,0,0,^8} 



{l,l,0,t9} 



{0,0,1,^10} 



{0,1,0,^11} 



{l,0,0,il2} 



{1,1,0,^13} 



{0,0,1,^14} 



{0,l,0,il5} 



TABLE II: Letter analysis of a command (gate) U 



I® I) 



2\/2 



(^r+af ) + 



~^2~y y y ' z <Jx ) 



4. A conditional Phase Shift. This gate provides a 
conditional phase shift c"^ on the second qubit 
(adds a phase to the second qubit): 



PhShift = 



10 
10 
10 

e*"^ 



-hi(3e-4 +e*$)I®I 
In our analysis we consider this phase to be 



We should notice that the succession of switch states 
follows a cyclic pattern. This regularity might facilitate 
the manipulation of the coupled junction device. 

Let us now give some numerical simulations of the 
proposed model for some fundamental quantum gates 
essential for the quantum computation (we use these 
gates multiplied with a proper constant because the cor- 
responding matrices should be elements of the SU{A) 
group) . These gates are : 

1. The CNOT gate. Probably the most important 
gate in quantum computation: 

10 
10 
1 
10 



CNOT 



(1) (2) , (1) 



J2) 



1(8)1) e SU{A) 



In the simulations of this paper, the energies are as- 
sumed to take the following values: 



^;c = 2.5 7^ = 3.4510"^^ J, 



El 



Er 



Q.IK = 0.138 10"^^ J, 



Q.IK = 0.13810"^^ J 



i.e. the time scale is of the order of 10^^^ sec . The nu- 
merical value corresponding to the idle state is chosen to 
conform to the available experimental data [^o| , and to 
keep in the range of different experimental propositions 
|]l5| , [1^ , |l^ |l^, |l9|. Each gate or command J/gate is ap- 
proximated by an evolution operator U{ti,t2, . . . ,^15) of 
the form (^2|). This simulation is equivalent to the anal- 
ysis of the command C/gatc to letters in conformity with 
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Table || The efficiency of our simulation is defined by a 
test function, ftest- It is a function of 15 time variables: 



ftest{tl,t2, ■ ■ ■ , tis) — 

= E l(f/gate)y-(C/(tl,t2, 



= u, 



gate 



-u\\^ 



Actually, ftest is the norm deviation of our simulation. 
The optimum, is obviously the nullification of this norm, 
ftest = 0. In fact we apply a minimization procedure and 
we calculate the time values, which minimize ftest- The 
numerical results are shown in Table [II . 



In our numerical examples we use an approximation of 
the time parameters to the fourth decimal digit. Respec- 
tively we calculate the value of the test function. Tak- 
ing into consideration three more decimal digits the test 
function ftest attains values of the order of 10~^°. It is 
a matter of intensity of the numerical algorithms used to 
find the minimum of the test function {ftest — 0) and 
convention of the number of decimal digits of the time 
parameters to succeed the optimal simulation . Indeed 
time parameters can not be determined with absolute 
precision in an implementation scheme for quantum com- 
putation. 

The construction of the one-qubit gates with the two- 
qubit Josephson device of Fig|| is possible. That is gates 
of the form I (g) and W (g) I, where W G SU{2), Fig |, 
which are simulated by the same device. 



W 



U 



W 



W(»I 



<w 



FIG. 3: Two-qubit and one-qubit gates in two qubit networks 



The construction scheme comes as follows: 



_ ^-itiHi -it-iHi -it-zHi -itiHi 



(14) 



where Hi and H/^ are special forms of the Hamiltonian 
(0), rewritten in the idle basis and ti, . . . ,1^ the time 
duration of each step. Obviously: 



1 p „(2) 1 p „(2; 



It can easily be shown that: 

U ^e- 



(13) Setting the total time 



E4fc7r 



i=l 



AE 



the previous relation is written as: 
U = 

The right hand side of the last relation is a 2 x 2 SU(2) 
matrix depending on 3 independent time parameters 
ti, t2, ^3, since the fourth time parameter t^ is specified 
from the demand that total time is assumed to be fixed. 
By an appropriate choice of these three time parame- 
ters any gate U of the above form can be constructed. 
So any command of the form U = I g) IF, which corre- 
sponds to an one qubit gate can be constructed by at 
most four steps. Therefore any command I g) IF can 
be analyzed in at most four letters. We simulate nu- 
merically the proposed model for the following one-qubit 
gates, I (g (NOT), I O (Had), I ® (VNOT) audi (g PhS: 



1(8) NOT 



10 
10 
1 
10 



I (8 (iax) e SU{A) 




^Vnot= ^ 



/ e *4 e*4 

e*f e~*^ 


V 











e* 4 



3PhS = 



/ 1 \ 

e"^ 

10 

\o e*'^ y 



cos — 1 
2 



t sm — (T, 
2 



The case of the gates of the form U = W (E) I can be 
similarly treated by using basic Hamiltonians Hi and H^ . 



Hi = [^Ee ai'^ - \Ej ai^)) + (\Ee of^ - \Ej <^ 



(2) 



where 



H3 ~ —\Ej(Tx {^EcCTz 



AE 
2 



,(2) 



1 P ^(2) 



o(2) 
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CNOT 




SWAP 


|1, 1, U, 102.7757) 


r 1 1 A TAA C O TOT 

|1, 1, (J, 700.5872) 


|U, U, 1, loo.oiyoj 


fA A 1 OAK OOAAT 

|0, 0, i, zOo.zoyUj 


f f\ 1 A AAA Ai^l 'tT 

|0, 1, 0, 909.9ol7| 


{0, 1, 0, 139.7456} 


ri A A lOAAC^A/lT 

|1, U, U, 13U.U5(J4j 


ri A A 1AA 0001T 

|1, 0, U, 199.3881 j 


fi 1 A OAA 'yOOAT 

|i, i, U, oUU. { 


fi 1 A 1 OA Aaaa^ 
|1, 1, U, lo0.4yDD| 


fA All AO ncyoT 
|U, (J, i, i4^.yo/oj 


fA A 1 1 1 K OA/l'yT 

|0, 0, 1, 110.894/ j 


TA 1 A 1A1 r\CO A~\ 

|0, 1, 0, 101.u5o4| 


TA 1 A IIAAC^O/lT 

{0, 1, 0, 110.9584} 


fi A A AAA CCiC\'\ T 

|1, U, 0, 90U.5o91| 


fi A A OAACttA/lT 

{1, 0, 0, 200.5504} 


{1, 1, 0, 151.529o| 


{1, 1, 0, 120.2794} 




/DDI 784 nnnsl 


{0, 1, 0, 161.0839} 


{0, 1, 0, 798.0702} 


{1,0,0,901.9591} 


{1,0,0, 129.1358} 


{1,1,0,699.3097} 


{1,1,0,501.6780} 


{0,0,1,191.4272} 


{0,0,1,130.0444} 


{0,1,0, 101.2086} 


{0,1,0,160.2219} 


/te.t = 2.9 X 10"" 


/test = 9.5 X 10"'' 



OFT, 




Phase Shift (0 = f ) 


|1, 1, 0, 710. 2581 j 




|1, 1, U, 110.71lDj 


fA A 1 AO/1 o/?oc:T 
|U, 0, 1, U84.8DO0 j 




fA A 1 AAO ^JOIOT 

|0, U, 1, 902.D81oj 


fA 1 A 1 ttA /I OATT 

{0, 1, 0, 159.4397} 




fA 1 A 1 OA OAOOT 

{0, 1, 0, 120.3082} 


ri A A 1 /I O C^^OAT 

|1, 0, 0, 142.5689) 




ri A A OAT OO/IAT 

|1, U, U, o97.824Uj 


fi 1 A 1 OO OT^iAT 

|1, 1, 0, loo. 2 (b\J\ 




fi 1 A 1 AA /lAAOT 

|1, 1, U, 109.4998} 


fA A 1 atio qkakT 
|U, 0, 1, DOO.OOOO j 




fA A 1 1 TK OIACT 

|0, U, 1, 1 (^o.819oj 


fA 1 A 100 /lOOCT 

{0, 1, 0, 133.4825} 




TA 1 A C01 ooaaT 

{0, 1, 0, 521.3209} 


ri A A AO /I OOOOT 

{1, 0, 0, 924.2883} 




fi A A lOOCACOT 

{1, 0, 0, 122.5053} 


{1, 1, 0, 173.8055} 




{1, 1, 0, 102.2305} 








{0, 1,0, 701.2262} 




{0, 1, 0, 108.9077} 


{1,0,0, 131.3048} 




{1,0,0, 198.7137} 


{1, 1,0,849.6562} 




{1,1,0,159.0513} 


{0,0,1,128.8483} 




{0,0,1,888.0009} 


{0,1,0, 150.0286} 




{0,1,0,100.7132} 


/test = 1.0 X 10" ' 




/test = 1.2 X 10"' 



TABLE III: Letter analysis of Two-Qubit Gates 



where 



IV. SUMMARY 



Ej 



The equivalent construction scheme is: 

jj _ ^—itiHi ^-itsHs Q~it2Hi ^—itiHi 

Obviously, 



(15) 



U 



^g— itfot Pz 



Setting the total time 



i=l 



Ak-K 



the previous relation is written as: 



U ■ 



-p.e-'*3— Tg- 



.1 



It is apparent that the numerical results for a simulation 
of an onc-qubit gate should be the same regardless of its 
form, 1®W or W So the numerical results concern- 



ing the time parameters presented in Table IV should be 



the same for the simulation of the corresponding gates 
NOT ® I, Had ® I, VNOT ® I, PhS ® I. 

Here it must be noticed that usually in order to 
achieve the construction of one qubit gates a more com- 
plicated technique is usually proposed. In this method 
the Josephson junctions are "neutralized" by appropri- 
ate annihilation of the tunneling amplitude Ej by using 
SQUID techniques P, p7f. 



The traditional approach to quantum computing is the 
construction of elementary one-qubit and two-qubit gates 
(universal set of quantum gates) which are connected by 
quantum connections and can represent any quantum al- 
gorithm pl[| . A different view is employed in the present 
paper, proposed in |^ under the name of encoded uni- 
versality. According to this, we do not force the system 
to act as a predetermined set of universal gates connected 
by quantum connections, but we exploit its intrinsic abil- 
ity to act as a quantum computer employing its natural 
available interaction. 

Thus, any one-qubit and two-qubit gate can be ex- 
pressed by two identical Josephson junctions coupled by a 
mutual inductor. This can be realized by a finite number 
of time steps evolving according to a restricted collection 
of basic Hamiltonians. These Hamiltonians are imple- 
mented using the above system of junctions by choosing 
suitably the control parameters, by switching on and off 
the bias voltages and the mutual inductor. The interac- 
tion times of the steps are calculated numerically. 

The values of the switches together with the values 
of the time steps may constitute the quantum machine 
language. Each command of the language consists of a 
series of letters and each letter of a binary part (the val- 
ues of the switch characterizing the Hamiltonian) and a 
numerical part (the interaction time). 

The generalization to A^-qubit gates is currently under 
investigation. In this case we need N + 2 basic Hamiltoni- 
ans in order to represent the corresponding A^-qubit gate. 
The structure of commands is an open problem. Each 
com man d can be obtained by 2^ — 1 letters see Proposi- 
tion A. 3 in Appendix ^ The mathematical foundation 
of this conjecture will be studied in another specialized 
paper. However by using the techniques described in 
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1(g) NOT 

{1,0,0, 133.6621} 
{1,1,0,104.2929} 
{1,0,0, 102.2461} 
{1,1,0,111.8267} 

ftest = 1.5 X llf^ 



NOT(g)I 
{0,1,0,133.6621} 
{1,1,0,104.2929} 
{0,1,0,102.2461} 
{1,1,0,111.8267} 

ftest = 1.5 X icr^ 



1(g) Had 

{1,0,0,113.3151}' 
{1,1,0,111.1253} 
{1,0,0,109.4076} 
{1,1,0,118.1799} 

ftest = 1.5 X 



hg)I 

{0,1,0, 113.3T5Ty 
{1,1,0,111.1253} 
{0,1,0,109.4076} 
{1,1,0,118.1799} 

ftest = 1.5 X 



I g) VNOT 
{1,0,0,115.7315}" 
{1,1,0,114.3695} 
{1,0,0,100.0236} 
{1,1,0,121.9033} 

ftest = 3.4 X llr^ 

VNOT g) I 

{0,1,0,115.7315}" 
{1,1,0,114.3695} 
{0,1,0, 100.0236} 
{1,1,0,121.9033} 

ftest ^ 3.4 X l"0^ 



Ig)PhS 

{1,0,0, 114.9978}" 
{1,1,0,109.2394} 
{1,0,0, 115.3838} 
{1,1,0,112.4068} 

ftest = 1.5 X 



PhSg)I 

{0,1,0,114.9978} 
{1,1,0,109.2394} 
{0,1,0,115.3838} 
{1,1,0,112.4068} 

ftest = 1.5X 



TABLE IV: Letter analysis of one-qubit gates in two qubit networks 



the number of letters can be reasonably reduced in the 
A^-qubit case. The application of the same methodology 
for other devices as quantum dots and NMR are under 
investigation. 



APPENDIX A: MINIMAL GENERATING SET OF 
THE SU(2'^) - ALGEBRA 

Let us consider the su(2) algebra in the adjoint repre- 
sentation. This algebra representation is a three dimen- 
sional vector space with basis the 2x2 Pauli matrices: 















(Jz = 













therefore 



and 



su(2) = span(cr^,(Tj,,(T^) 



i 
-i 



This su(2) algebra in the adjoint representation is gen- 
erated by the following 2x2 hermitian matrices: 



1 
1 



and (Tt 



1 

1 



because 



[<7z,crx 



2icr„ 



The adjoint representation of the algebra su(2^) is the 
vector space spanned by the 15 matrices 

su(22) = span(af \ , , t,j - 1,2,3) (Al) 

where 

o-f ) = tJi (g) I, fjf ^ =I(E)a„ (jf ^ = fJi (g) a J 

the adjoint representation of the su(2^) algebra can be 
generated by linear combinations and successive commu- 
tations of the following 4 elements: 



This is indeed true because all the elements of the basis 



(Al) can be generated by repeated commutations of the 
elements (|A^), because for fc = 1, 2, 3: 



.(1) 



^(fe) ^(fe) 

(Jz , Oy 



(A3) 



Oz , 



.(2) 



The elements cr^^V^^^ can be generated by commutating 
the generators am with a'\f \ One illustrative exam- 
ple is the const ruct ion of the element ai^Vi^"*, by using 
the generators ( A2): 



(J X (J z 



16 



(J X , 



(1) (1) (2) 



^ Z 1 



.(2) 



(1) (1) (2) 



(A4) 

In the case of the adjoint representation of the algebra 
su(2'^), we can work following a similar methodology. The 
adjoint representation algebra su(2^) is a vector space 
spanned by the following 63 matrices: 



1(23) 



enanfrr^^^ ^(2) _(3) ^(1)^(2) 
span^cTj , cr- , (7^ , cr- , 

afVf af\z,j,fc = l,2,3) 



(A5) 



where 
^(1) - 



1(g) I, 



.(2) 



I (g CTi g) I a. 



(3) 



Ig)I( 



The above elements can be generated by repeated com- 
mutations of the following 5 matrices: 



.(1) 



Oz , 



,(2) 



Oz , 



"v "y 



cri^^ , and 



(A6) 



The linear terms al can be easily generated by formu- 

Jk) (e) 

3 



a'f^ , af'^ , cr(^) -I- af^ , and a^y'^ af^ (A2) las as in equation (Q. The quadratic terms af' ay are 
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generated by manipulations slightly more complicated 
than in the case of equation (A4). Let us take the ex- 



ample of the generation of the element cy^^af^^ then we 
must perform the following commutation actions: 



,(1)^(2) 



(2)^(3) _ 



da; CT 



(2) mm 



y 



.(1)^(3) 



.(2)^(3) 



U X Ox 

„(2) _ i 



(1) (1) (2) 



O z 1 Ox 



,(2) 



(A7) 



.(2) ^(1)^(2) 



Therefore all the quadratic terms can be generated by 
the elements (A6). Let us now generate a cubic term 

of the algebra as the element u^z^ a'^jp This element 
is generated by the commutation elements and 
o'^'^ (J^\ which are generated previously: 



These relations prove that the su(4) algebra can be gen- 
erated by combinations and successive commutations of 
the four elements {Hi, H2, H3, H4}. Starting from this 
fact we can construct all the elements of the form: 

_Q—iH^t^^—iH^t2Q—iH\t\ 

(All) 

From the Baker-Campbell-Hausdorff (BCH) formula |2|, 
Sec. 2.15]: 

where 

h{A,B) = A + B + ^[A,B]- 

-^{[[A,B].B]-[[A,BIA])- 
-i-,{[BM. [A.B]]]-[A, [B,[B,A]]]) + ... 



one can calculate U in (All) by successive applications 
of the the above BCH formula starting from the left to 
the right, i.e. 



By induction we can prove the following proposition: 

Proposition A.l The adjoint hermitian representation 
of the algebra su{2^), i.e the set of hermitian traceless 
2^ X 2^ can generated by the algebra of Lie-polynomials 
of the set: 



N 



N 



A 



N 



V V 



(A8) 



The set of the generators has -I- 1 elements, we 
should point out that this number is very small than the 
number 4^ — 1, which is the dimension of the algebra 
su(2^). Therefore, large Lie algebras can be generated 
by using a relatively small number of elements. 

Let us now construct the group SU(2^). For the sake 
of simplicity we start the discussion with the SU(4) case, 
i.e with the set of unitary 4x4 matrices with determinant 
equal to 1. 

Let us consider four linearly independent elements, 
which are given by the formulas: 



Hi = i^c {<tP + 
H2 = -\Ej iai'^ 

rj _ 1 77. ^(2) 
Hi 



J2)^ 



(1) 



,(2)n 



2'^j{^x + CTx j 



\Ej {a. 

Ip ^(1)^(2) 



(2)n 



(A9) 



kEj (4 



.(2) 



Starting from this system we can reconstruct the ele- 



ments 



because: 



^1 ^ 2(gi-g3) 
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2 _ 2(Hi-H3-Hi) 
Ej 

2{H3+Hi-Hi-H2) 
El 



(AlO) 



U = exp[— TO15] where ui = ii-ffi, 

U2 = h{t2H2,Ui), 
U3 = h{t3H3,U2) 



Ul5 = h{ti5H3,Ui4) 

The elements Uk are complicated combinations of the el- 
ements Hi, H2, H^, H4, bracketed inside commutators, 
i.e are Lie-polynomials of the free associative algebra 
C{Hi, H2, H3, H4}. By definition the combinations and 
all the successive commutators of these elements gener- 
ate the algebra su(4), which has as a linear basis the 15 
elements given by equation (Al). Thus the general form 
of U15 is given by: 



Ul5 



i=l 



f(2)^(2) 



where the 15 coefficients 



Ji — Ji 1^1; 12,. 

and 

„(1,2) _ .(1,2). 
Ji - Ji 111,12, 



,tl5) 



(A12) 



are complicated functions of the parameters 
ti,t2, ■ ■ ■ ,ti5. Any element of the algebra su(4) is 
written as a linear combination of the elements of the 



basis (Al) and from the known functions (AIS) we can 
find the values of the finite time series ti, i2,-..,^i5- 
Then we have proved the following Proposition 

Proposition A. 2 The group SU{4) is given by the ele- 
ments of the form: 



U 



,Q-iHit5Q-iHiti^-iH3t3^-iH2t2Q-iHiti 



where Hi, H2, H3, H4 are some special elements of the 
algebra sm(4). The combinations and the successive com- 
mutations of these elements generate sw(4). 
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This proposition is a special form of the bang-bang con- 
troUabiUty theory for SU(4) matrices |^. The above de- 
composition of the SU(4) matrices is a generahzation of 
the Euler decomposition of the SU(2) matrices. In ||] 
another decomposition is proposed based on the Cartan 
decomposition of SU(2") matrices. In the Cartan de- 
composition a choice of orthogonal basic Hamiltonians is 
used. In the present decomposition we do not consider 
an orthogonal set of basic Hamiltonians but the form 
of the Hamiltonians is imposed by the physical system 
under consideration, fulfilling the conditions of encoded 
universality |j 

Proposition A. 3 Let Hi, H2, H^, H4, are elements 
of the algebra su{2^) such that the successive commuta- 
tions of these elements generate su{2^). Then any el- 
ement U belonging to the group SU{4:) is given by the 
relation: 



A similar Proposition can be stated in the case of the 
general problem related to the SU(2") group. The above 
Proposition concerns the problem of controllability of 
spin systems, in the context of the Cartan decomposi- 
tion technique. This problem was solved |^ and other 
studies of the same problem by different techniques have 
been recently proposed |23, The general problem 



can be formulated in a different way. From the theory 
of universal gates ||l| and the papers on the control of 
the molecular systems ^ it is well known that the 
811(2") can be decomposed into simpler matrix factors 
with SU(2) and SU(4) structure. That means that the 
one and two qubit gates are universal ones. A systematic 
study of this technique is under current investigation. 
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